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We use magnetic force microscopy (MFM) to measure the local penetration depth A in 
Ba(Feo.95Coo.o5)2As2 single crystals and use scanning SQUID susceptometry to measure its tem- 
perature variation down to 0.4 K. We observe that superfluid density ps over the full temperature 
range is well described by a clean two-band fully gapped model. We demonstrate that MFM can 
measure the important and hard-to-determine absolute value of A, as well as obtain its temperature 
dependence and spatial homogeneity. We find ps to be uniform on the submicron scale despite the 
highly disordered vortex pinning. 



The magnetic penetration depth A, one of the two fun- 
damental length scales in superconductors [l] , character- 
izes many fundamental properties. It evaluates the phase 
stiffness of the superconducting state by the temperature 
rpmax ^ 1/A^ at which phase order would disappear [2 . 
It also determines the superfluid density ps = 1/A^, the 
number of electrons in the superconducting phase. How- 
ever, its absolute value is notoriously difficult to measure, 
especially in samples that may have either intrinsic or ex- 
trinsic inhomogeneity. In this letter, we will report a new 
technique to measure A by magnetic force microscopy 
(MFM). The advantage of using local probes over bulk 
techniques is that it allows us to study the sample homo- 
geneity. We implement this technique to determine ps in 
a iron-pnictide superconductor Ba(Feo.95Coo.o5)2As2. 

Iron-pnictides superconductors have been under exten- 
sive study since their recent discovery [3 . The high tran- 
sition temperature [4 , the proximity to a magnetic state 
and the existence of multiple conducting bands 
m |9] combined to make it difficult and interesting to 
resolve key issues like the superconducting order param- 
eter (OP) symmetry flOjIII], the pairing mechanism [12 
and the role of impurities and inhomogeneity [13] . Those 
problems can be studied by measuring pg. When the gap 
has nodes, ps{T) varies as a power law in T at low T, 
as demonstrated in YB2Cu307_5 [HlIIS], while a fully 
gapped OP gives a low-T exponential dependence [16]. 
Since it is difficult to determine A, its temperature vari- 
ation AA(T) = A(T) — A(0) is often measured, which 
follows the same temperature dependence as ps at low 
T. Sometimes this approach is sufficient, e.g. linear AA 
in clean LaFePO over a wide temperature range provides 
strong evidence of well formed line nodes [TTl [18] . How- 
ever, in the Ba-122 family, a steep power-law A A was 
obtained in the Co-doped compounds [I9l [20] while an 
exponential ps was measured in the K-doped materials 
[2T] . The question waiting for clarification is whether 
different dopants lead to different OP structure. AA mea- 
surement can not infer OP symmetry except for T <C Tc, 



but for multi-band pnictides, the low-T regime may be 
dominated by the small-gap regions of the Fermi surface 
and may be altered by interband impurity scattering [22] . 
It is thus important to measure the absolute value of A 
to determine ps over the full temperature range. 

In this paper, we measure the local AXab{T) 
and Xab{T)^ the penetration depth for screening cur- 
rents flowing in the a-b planes, in electron-doped 
Ba(Fei_a;Coa;)2As2 single crystals {x ^ 0.05, Tc = 18.5 
K, grown from self- flux [6 ) from T = 5 K to Tc by 
magnetic force microscopy (MFM) [Fig.[l]. We also use 
scanning SQUID susceptometry (SSS) [23 to measure 
AXab{T) down to 0.4 K. We find that ps can be well 
described by a two-band fully gapped OP over the full 
temperature range. We also use MFM to image and ma- 
nipulate vortices to measure the homogeneity of Xab{T) 
and the flux pinning force. We find that ps is uniform to 
within 10% or better, although vortex pinning is highly 
inhomogeneous . 

In our MFM, a sharp magnetic tip at the end of a flexi- 
ble cantilever faces the crystal surface, which is parallel to 
the a-b plane. By measuring the shift in the cantilever's 
resonant frequency [24 , we determine dF^/dz [25 , where 
F is the force between the tip and the sample, and z is 
along the tip magnetization direction and is normal to the 
cantilever and to the crystal a-b surface. dFz/dz changes 
abruptly within a few nanometers of the surface, allow- 
ing precise determination of the tip-sample separation 
z. In the Meissner state, the tip-superconductor inter- 
action can be approximated by the magnetic interaction 
between the tip and its image mirrored through a plane 
at z = —Xab (Fig. [1^ inset) [26]. This local levitation 
force is determined uniquely by z -\- Xab{T) for z ^ A^^ 
(Ac does not enter for any source field above a smooth, 
infinite ab surface) [27 . Thus, changing T at constant z 
offsets a dFz/dz curve along the z-axis by AXab{T). To 
acquire the data labeled as MFM AA in Fig. [2| we park 
the tip at z = 500 nm, change T and acquire dFz/dz. 
The z offset required to match dFziT) / dz with a refer- 
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ence curve at T = 5 K gives \ah{T) — Aa5(5K). Using a 
similar method for data acquired by SSS in a ^He refrig- 
erator [18], we extend measurements of AA down to 0.4 
K on two nominally identical samples. The SSS results 
match the MFM results over the common temperature 
range. By using local scanning probes, we reduce the 
influence of the complex topography around the sample 
edges [27 . 

Figure [2] shows that AXab{T) increases very slowly with 
T at low T, inconsistent with the the linear dependence 
that would be expected for line-nodes. The same behav- 
ior appears at three different locations on two samples 
with SSS and at four different locations with MFM on a 
third sample. Between T = 0.02Tc and OATc AXab{T) 
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varies by about an order of magnitude less than has been 
reported for a similar sample using a bulk technique [20] . 
At low T, AXab{T) can be described by either a two- 
band fully-gapped model or by a power law with a small 
coefficient as described below. 

We also extract Xab{T) by modeling the tip- 
superconductor interaction, with the magnetic tip as 
a sharp, single domain cone, truncated a at distance 
ho = 400 ± 20 nm from its apex as shown in Fig. [l^. 
Within the model, the z-dependence of dF^/dz is given 
by: 



dF,{z,T)/dz 
. / 1 



dF,{z,T)/dz\ 
ho 



(1) 



Z + Xab{T) {z + Xab{T)y 2 {Z + Xab{T)y 



where A is determined by the tip shape and the coat- 
ing. The value A = 78pN from fitting at T < Tc/2 is 
consistent to within 30% with the magnetic moment ex- 
pected from the nominal iron coating on the tip, and 
with that inferred from the tip- vortex interaction [25] . 
We record dFz/dz as a function of z and T and extract 
Xab at many temperatures by fitting to Eq. [ijwith A and 
ho fixed and Xab and dFz/dz{oo^T) allowed to vary sep- 
arately for each T. The fit works well for all T (Fig. \\\. 
The resulting values of Xab{T) are shown in Figure [2| 
with label "MFM A" and agrees well with the model- 
independent AA. If we consider only statistical errors, we 
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FIG. 1: Technique to measure A and A A by MFM from Meiss- 
ner repulsion (a) and vortex imaging (b,c). (a) z dependence 
of dF^/dz (blue symbols) at T = 5, 12 and 18 K and the 
fit to the truncated cone model (red dashed line). (Inset): 
Sketch to illustrate that the tip-superconductor interaction in 
the Meissner state can be approximated by the interaction be- 
tween the tip and its image mirrored through a plane (dashed 
line) Xab below the surface of the superconductor (solid line) 
when z ^ Xab- Comparing the curves provides AXab inde- 
pendently of the tip model. Fits give Xab{T) at T = 5, 12, 18 
K to be 0.33, 0.37, 1.10 /am. (b, c) Images of two vortices 
{z = 400 nm) at 5 K (b) and 10 K (c). The shapes and am- 
plitudes depend on both the magnetic field from vortices and 
the tip structure, but the similarity shows that both the spa- 
tial variation and the temperature-induced change of Xab are 
small. (cl,e) Scanning electron microscopy images of the tip 
before (d) and after (e) the measurements. Also shown are the 
truncation distance ho — 300 ±30 nm in (d) and /lo = 400 ±20 
in (e). An accidental crash during the measurement changes 
the truncation distance ho from 300 ± 30 nm (d) to 400 ± 20 
nm (e). Despite the crash, dFz/dz curves taken before and 
after the crash give the same Aa6(5K) to within 10 nm. 
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FIG. 2: Normalized superfiuid density ps{T)/ps{0) = 
Xab{Of/Xab{Ty vs. T. We determine AAa6(T) by MFM 
(squares) and by SSS (diamonds) from measuring the change 
in the diamagnetic response at fixed height. These values 
are offset to match the absolute value of Xab{T) obtained by 
fitting the MFM data to the truncated cone model (circles). 
The green solid line shows a fit of the two-band s-wave model 
discussed in the main text (Ai = 2.6Tc, A2 = 0.8Tc, x = 0.88 
and a = 1.4). The width of the dashed band reflects the uncer- 
tainty in Xab{0). Inset: AAab vs. T at low T. Black dashed 
line: one-gap s-wave model with a = 1.5 and Ao = 1.95Tc. 
Magenta dashed line: AXab{T) = cT^"^ ( c = 0.14nm/K^-^). 
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obtain Aa6(5K) = 325 ± 5 nm with 70% confidence inter- 
val. However, the systematic error from the finite width 
corrections of the tip-geometry is 5%. In addition, the 
±20 nm uncertainty on ho leads to 74 pN < A < 81 pN by 
bootstrapping. The extremals of A and ho gives ±35 nm 
systematic error on Aa6(5K). Thus, adding the two main 
sources of systematic error, we find Aa6(5K) = 325 ± 50 
nm. 

Knowing Xab{T) gives ps over the full temperature 
range (Fig. |2|. The fact that ps does not saturate 
at low T is inconsistent with a single-band isotropic 
gap. A two-band fully gapped OP, which was pro- 
posed theoretically |TT] and tested experimentally 
O m], describes the data weh (Fig. [2|. In the model, 
Ps{T) = xpi{T) ± (1 - x)p2{T): pi,2(T) are the super- 
fiuid densities in bands z = 1, 2, with gaps A^(T) = 

A^(0) tanh ^ a ^o) \/^^ ( ^ ~ ' describes the rate of 

Ai{T) increasing upon cooling from Tc [IB^. Our fit (tak- 
ing into account the systematic error on Aa6(5K)) gives 
Ai(0) = 2.5±0.3Tc, A2(0) = 0.70±0.1Tc, x = 0.89±0.06 
and ai = 1.45 ± 0.4 with a2 = 1. The value of ai sug- 
gests that pairing is likely to be more complicated than 
phonon- mediated weak coupling [TJl |28], which would 
give a = 1. The magnitude of Ai^2(0) is consistent with 
the scaled down values deduced from optical spectroscopy 
on similar materials with higher Tc [H |9j . At low T a 
power law cT^ where n = 2.2 and c = 0.14 nm/K^*^ 
also fits the data. The dominant sources of errors are the 
calibration accuracy of the scanner, thermal drift, and 
the breakdown of the assumption of z ^ Xah-, which to- 
gether would bound c between 0.12 and 0.18 nm/K^-^. 
The small coefficient is inconsistent with that previously 
reported [19]. We rule out a nodal OP model since the 
impurity scattering rate required for such a model [29 to 
match our data is much higher than that reported in pre- 
vious works on d-wave cuprates with deliberately added 
impurities ^^|3T]. Instead, we interpret this weakened 
exponential behaviour of AA^b and ps{T) from 0.4 K all 
the way to Tc as strong evidence for two full gaps, con- 
sistent with the extended s-wave OP [TOl [32] . 

We repeated the touchdown measurement at four posi- 
tions separated by around 10 pia and obtained Xab{T = 
5 K) = 325 nm, 330 nm, 325 nm and 330 nm. This result 
suggests that Xab is uniform across the sample. 

A second test of uniformity is afforded by measuring 
the local Tc by mapping the lowest T at which we can- 
not detect Meissner levitation by MFM (sensitivity cor- 
responds to Xab{T) > 3 pm) or diamagnetic response by 
SSS (sensitivity corresponds to Xab{T) > 20 pm [33^). We 
find the variation of Tc to be less than 0.5 K throughout 
the range of 10 x lOpm^ by MFM and 200 x 200/im2 by 
SSS. 

Vortex imaging provides a third test of ps{T) unifor- 
mity. To this end, we cool the sample in an external 
magnetic field and scan the tip at a constant height z 



above the surface at 5 K. All vortices appear very simi- 
lar (Fig. |3^), indicating that the spatial variation of Xab 
is limited. The convolution of the tip and the vortex 
field makes it difficult to extract A^^ from the vortex 
imaging. Instead, we calculate the normalized curva- 
ture at each vortex peak to quantify the spatial variation: 

C = md.x{dF,/dzy^det ( ^'£fg^f ^ ) {h j run over 1, 2 

and xi = X, X2 = y)' The length-scale, C~^/^, charac- 
terizes the spatial extent of the magnetic field from each 
vortex (Fig. [3^). The scatter (Fig. [sJd) of the normal- 
ized C~^^^ (±8%) at constant z gives an estimate for the 
spatial variation of A^^. 

In contrast to the uniform p^, vortex pinning is very 
inhomogeneous. Vortices do not form an ordered lat- 
tice when field-cooled in fields up to 13 mT, the highest 
field that allows us to resolve individual vortices in this 
material. Instead, vortices always appear in the same 
regions when we thermal cycle in different fields using 
different cooling rates (Fig. |4^). This behavior suggests 
inhomogeneous pinning. To measure the pinning force 
distribution, we use the MFM tip to drag individual vor- 
tices and to convert the recorded dFz/dz to the required 
force [25]. We measure two different forces (Fig.|4]3): the 
force for dragging the most weakly pinned vortex, Tmin, 
a measure of the smallest pinning force (Fig. |4]:); and 
the force for dragging all of the vortices (usually ^10) 
in a field of view, Ttyp, a measures of the typical pin- 
ning force (Fig. [iji). In this sample 2 < Ttyp/Tmin ^ 4. 
Ttyp ~ 18 pN at 5K, corresponding to a critical current of 
Jc ^ 80 kA/cm (Tc = Jc^o^i where is the fiux quan- 
tum, d = 10pm is the sample thickness), consistent with 
the value from bulk measurement of an optimally doped 
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FIG. 3: Spatial uniformity of Xab from vortex imaging at 
5 K. (a) Image of vortices at T = 5K, z = 125 nm and a 
vortex density giving 3.5 mT. (b) The normalized length-scale 
associated with each vortex peak in (a): C~^^^ j < C~^^^ > 
(<> denotes the mean). We find C"^/^ = 250(1 ± 0.08) nm 
(70% confidence interval). 
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doped iron-pnictides may be suppressed by the competi- 
tion with non-superconducting phases. MFM allows us 
to obtain the superfluid density and to map its spatial 
variation down to the submicron scale. This capability 
may be useful to study how different phases compete for 
charge carriers. 
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FIG. 4: Inhomogeneous vortex pinning, (a) Image of vor- 
tices at T = 5 K, z = 80 nm and B — 9.5 mT, overlaid by 
the vortex positions (dots) in Fig. [3^ and the boundary of 
that scan (black frame). The vortex configuration is highly 
disordered. Vortices avoid the same regions in both scans, 
taken days apart and many thermal cycles apart, (b) Local 
critical current (left ordinate) and the depinning force (right 
ordinate) vs. T. The comparison of minimum and typical 
values implies inhomogeneous pinning, (c) Image of vortices 
atT = 5K, z = 120 nm showing that Fmin only moves 
the vortex at the bottom, (d) Image of moving vortices at 
T = 14. 5K, z = 430nm showing that Ftyp allows us to drag 
all vortices a distance of several microns. 

sample [34^- Even at i^typ? vortices do not follow the 
tip all the way, indicating the existence of pinning forces 
larger than Ftyp- In fact, Ftyp is still at least an order 
of magnitude smaller than the force required to stabilize 
vortices in the dense clusters we see (the vortex-vortex 
interaction for a pair separated by 400 nm corresponds 
to a current density of 3 MA/cm^). We do not detect 
any correlation between pinning and superfluid density, 
suggesting that strong pinning exists without affecting 
superconductivity on the scale of Xab- The ability to mea- 
sure the absolute value of the penetration depth despite 
a disordered vortex configuration is important, since the 
most commonly used method, muon-spin-rotation [35] . 
assumes an ordered vortex configuration. 

To conclude, by measuring Xab{T) and AXab{T) locally 
we find that underdoped Ba(Fei_a;Coa;)2As2 {x ~ 0.05) 
has homogenous ps whose temperature dependence can 
be described by a two-band fully-gapped OP. This result 
provides thermodynamic evidence for fully gapped mod- 
els such as the proposed extended s-wave model [lOl [TT] 
for Co-doped 122 pnictides and shows that it has the 
similar OP structure as the K-doped, despite the differ- 
ent dopants and substitution cites. We obtain Xab{^) = 
325 ± 50 nm, which gives TJ^^^ = A{hc)'^a/{167re'^ X'^) ^ 
260 K, where a = \/7r^c5 = 1-1 nm ^ and A=2.2 in 
the three-dimension limit [l [34]. T^^^^ > hinting 
that phase fluctuations are not as important here as in 
the underdoped cuprates [2|. Instead, Tc in the under- 
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